Colloidal brazil nut effect in sediments of binary charged suspensions 
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Equilibrium sedimentation density profiles of charged binary colloidal suspensions are calculated 
by computer simulations and density functional theory. For deionized samples, we predict a colloidal 
"brazil nut" effect: heavy colloidal particles sediment on top of the lighter ones provided that 
their mass per charge is smaller than that of the lighter ones. This effect is verifiable in settling 
experiments. 
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Binary systems of granular matter separate upon shak- 
ing in gravity, so that the larger particles lie on top of the 
smaller ones even if they are heavier and denser than the 
latter. This is due to a sifting mechanism in which tiny 
grains filter through the interstices betv^^een the large par- 
ticles which is well-known as "brazil nut" effect: in a jar 
of mixed nuts or in a package of cereal, the largest species 
rises to the top P, 01 • This clearly distinguishes granular 
matter from ordinary fluids where the rising species is 
controlled by Archimedes' law. Understanding the full 
details of the brazil nut effect is still a problem; recently 
even a reverse brazil nut effect of large light grains sink- 
ing in a granular bed has been predicted |^ |^ and 
verified in experiments j^. 

In this letter, we report on equilibrium density profiles 
of binary charged colloidal fluids ( "macroions" ) under 
gravity. Using extensive Monte-Carlo computer simula- 
tions of the "primitive" model of strongly asymmetric 
electrolytes and a density functional theory, we predict 
that the heavier particles sediment on top of the lighter 
ones provided the charge per mass of the heavier parti- 
cles is higher. In analogy to granular matter, we call this 
counter-intuitive phenomenon a colloidal brazil nut ef- 
fect. It is generated by the entropy of the microscopic 
counterions in the solution, which are coupled to the 
macroions by strong Coulomb binding. Clearly, though 
this effect is qualitatively similar to the granular brazil 
nut effect insofar as heavy particles are on top of lighter 
ones, its physical origin is different: first, the particle 
charge (and not the size) is crucial. Second, the col- 
loidal brazil nut effect is a pure equilibrium phenomenon 
while the granular brazil nut effect happens intrinsically 
in non-equilibrium. The colloidal brazil nut effect can be 
verified, e.g., in depolarized-light scattering or real-space 
experiments on sediments of strongly deionized binary 
charged suspensions H, . Similar techniques have been 
used to measure one-component colloidal density profiles 
Il0| where deviations from the ideal barometric law 
are still an ongoing debate 12, 13, 14, 15, 16] . 

We simulate the asymmetric "primitive" model of bi- 
nary charged suspensions in which the solvent only enters 
via a continuous dielectric background with permittiv- 
ity e but all charged particles (two species of negatively 
charged macroions and microscopic counter- and colons) 



are treated explicitly at constant temperature T. With 
Zie, Zie, —qe and txi, CT2, (Tc denoting the charges and 
the diameters of the two colloidal species and the mi- 
croions, the interaction between the charged particles is 
given as a combination of Coulomb forces and excluded 
volume of the hard particle cores. Here, e is the electron 
charge. For simplicity we assume that the co- and coun- 
terions of the salt solution have the same valency and the 
same hard core diameter CTc. We consider a finite system 
of Ni , N2 charged macroions and corresponding number 
of counterions (fixed by global charge neutrality) plus 
salt ions of bulk concentration . The simulation box is 
rectangular with lengths = Ly and Lz = 32 L^, in the 
three different spatial directions with periodic boundary 
conditions in x and y direction and finite length in z- 
direction. Hard walls are placed at z = 0, and gravity 
with acceleration g points along the —z direction. Only 
the colloidal particles with their buoyant masses mi and 
TO2 are subject to gravity, whereas the microions are not. 

A Monte-Carlo simulation is performed in the canon- 
ical ensemble with the long-ranged Coulomb interaction 
treated via Lekner sums [l3- Typically lO'^ Monte Carlo 
moves per particle were performed for equilibration and 
it took an additional 10^ Monte Carlo moves per particle 
to gather statistics. Finite system size effects were care- 
fully checked by changing all lateral linear dimensions by 
a factor of 4. This means that we have changed the total 
number of colloidal particles A^i -I- N2 in the range of 12- 
200. We have calculated the inhomogeneous z-dependent 
averaged density profiles pi{z), P2[z), ps^z) and P4(z) of 
the two macroions and the counter- and colons. Data are 
shown for the largest system size where iVi ^ N2 — 100. 

Besides the ratios Zi/q, Z2/q, mi/ 1712, <7i/o'2, tTc/0'2, 
the system is characterized by two partial area den- 
sities Ui = Ni/LxLy {i = 1,2), the Bjerrum length 
\b = q^e^ /ekT [kT denoting the thermal energy), the 
gravitational length I2 — kT/m2g of the second parti- 
cle species, and the salt concentration. In order to re- 
duce the parameter space, we have assumed throughout 
the simulations monovalent microions {q ~ 1), and the 
same hard core diameter of the macroions a — ai — 0^2 
which serves as a natural length scale. We further fixed 
Z2 = 15, As = (Tc = c^/128 and m = n2 = 0.1/a^. This 
corresponds to typical parameters for low-charge aque- 
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FIG. 1: Density profiles of macro- and counterions for Z\ — 
45,30,25 (solid lines; top to bottom). The second colloidal 
component is shown as a dashed line. Counterion densities 
(thin lines) have been divided by Zi + Z2- For clarity, curves 
pertaining to different simulation runs have been shifted by 
5 X 10~* with respect to each other. The parameters are: 
m\/m2 — 1.5, ^2/0" ~ 10, Cs — 0. Inset: Semi-logarithmic 
plots of the same colloidal density profiles with the slopes 
predicted from theory. 



ous suspensions. We varied the colloidal charge Zi, the 
colloidal mass ratio mi/m2 (with mi > m2), the gravi- 
tational length £2, and the salt concentration. 

In the salt-free case, density profiles for the two 
macroions and the counterions are shown in Figure ^ 
for mi/m2 — 1.5 and three different macroion charges 
Zi = 45, 30, 25. For large heights z, the heavy particles 
(solid curve) are on top of the lighter ones (dashed curve) . 
This colloidal brazil nut effect is getting stronger for in- 
creasing charge asymmetry between the colloidal parti- 
cles: for the highest charge Zi = 45, the two macroion 
species are almost completely separated. To quantify the 
brazil nut effect, we define - in analogy to the ordinary 
brazil nut problem - a mean (or sedimentation) height 
of the two profiles via 



j'^ zp,{z)dz 
Pi{z)dz 



z = 1,2 



(1) 



By definition, the colloidal brazil nut effect occurs if hi > 
/i2, but there is no brazil nut effect (or a "reverse" brazil 
nut effect) in the opposite case, /ii < /i2. 

Let us now describe the properties of the density pro- 
files for increasing z qualitatively and put forward a sim- 
ple theory to describe the basic features. Close to the 
hard container bottom at z = 0, small correlation effects 
are visible as a density shoulder of the low-charge par- 
ticles while high charge particles are depleted from the 
wall. The reason is a combination of a pure interface ef- 
fect and Archimedes' principle. First, we have checked 
by simulation that depletion of high-charge particle per- 
sists for zero gravity. Second, by crudely mapping the 
interacting binary colloidal mixture onto one with effec- 



tive hard interaction cores, the high-charge particles will 
have a larger diameter and hence Archimedes' principle 
will lift the high-charge particles in the sea of small ones 
provided their mass density is smaller Il8'|. For our pa- 
rameters, however, the latter effect is small and confined 
to regions close to the wall. More importantly, as re- 
vealed by the semi-logarithmic plot in the inset of Fig- 
ure n there is an exponential decay of the two colloidal 
profiles for intermediate heights associated with two dif- 
ferent decay lengths for the two colloidal species. The 
decay length for the heavier but high-charge particles is 
larger than that for the low-charge particles. This gives 
the most significant contribution in the integral (1) of 
the mean height. Finally, for large z, there is cross-over 
towards another exponential decay involving the gravita- 
tional length £1 = kT/niig for the high-charge particles. 
It is important to note that the density profiles fulfill local 
charge neutrality throughout the whole sample except at 
the container bottom (z = 0) and at very large heights z. 
But even there the excess charge separated is small [Tq). 

Our theoretical explanation for the colloidal brazil nut 
effect is based on a simple density functional approach. 
The free energy J- per unit area, which is a functional 
of the inhomogeneous density fields pi{z) , p2{z) , p3{z) 
of the macro- and counterions, splits into the gravi- 
tational energy, the entropy of the three species and 
all Coulomb contributions. The latter are approxi- 
mated within a mean-field-type Poisson-Boltzmann the- 
ory El El El [Hill. Hence: 
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Here, A,y, (1/ = 1,2,3), are Lagrange multipliers which 
ensure that the overall colloidal densities per unit area 
equal the prescribed number densities, i.e. Jp°° dz pi{z) = 
Hi {i = 1,2). A3 is fixed by global charge neutrahty, 
and pt{z) — Ziepi{z) + Z2ep2{z) — qep^{z) is the total 
local charge density of the system. The functional J- is 
minimal for the physically realized equilibrium density 
profiles. 

We discuss two different cases of weak and strong 
Coulomb coupling subsequently. For weak Coulomb cou- 
pling, which is realized for very large heights z where 
the densities are extremely small, one may neglect the 
third term on the right hand side of Eq. (2). Then, the 
minimization of J- yields colloidal density profiles which 
follow the traditional barometric law pi{z) cx exp(— z/£i) 
{i — 1,2). Strong Coulomb coupling, on the other hand, 
will impose local charge neutrality pt{z) = which im- 
plies that the counterion density field is enslaved to that 
of the macroions. Minimizing T in this limit with respect 
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FIG. 2: Crossover towards the brazil nut effect. Shown are 
different charge asymmetries {Zi + 1)/{Z2 + 1) and mass ratios 
mi/m2 for which the brazil nut effect was (circles) or was not 
observed (crosses). The straight line shows the theoretical 
prediction. The parameter combinations were obtained by 
changing £2 and Z\ is the range of 4...20(t and 16 . . . 45, 
respectively. 



to the two colloidal densities only, again yields an expo- 
nential decay Pi{z) oc exp(— 7^2) (i = 1,2) but with in- 
verse decay lengths 7^ which are smaller than 1/ These 
decay lengths 7i turn out to be as follows: Let a be the 
index i for which the mass per charge ratio z"/^q+i 
minimal and let (3 be the index complementary to a, i.e. 
/3 = 1 if a = 2 and /3 = 2 if a = 1. Then 



and 
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(3) 



(4) 



with 7/3 > i;x(z^/q+i) — 7"- "^^^^ second case of strong 
Coulomb coupling will be realized for heights z where cor- 
relations between the ions are small (which justifies the 
mean-field approximation) but where local charge neu- 
trality is still valid. The physical reason for the much 
slower decay of the colloidal density profiles in the sec- 
ond regime results from the counterion entropy which 
tends to delocalize the counterions. However, since the 
macroions are coupled to the counterions due to the con- 
straint of local charge neutrality, they are lifted upwards 
together with their counterions. Assuming that the main 
contribution in the integrand of the right- hand-side of (1) 
comes from the second regime, the mean heights are given 
by ha = l/"fa E^nd hf} = I/7/3. A brazil nut effect occurs 
when a = 1 and P = 2. Hence the transition towards the 
brazil nut effect happens at zjq+i ~ zjjq+i ■ 

Let us now test the prediction of the theory against our 
simulation data. First, the slopes in the inset of Figure^] 
confirm the inverse decay lengths 7q and 7^ perfectly 



FIG. 3: Mean heights hi of heavy {+) and light (0) particles 
as a function of Zi for the parameters of Fig. The lines 
show the theoretical predictions. 



for a large range of intermediate heights; the theoreti- 
cal predictions for the slopes as given by Eqs. (3) and 
(4) are shown as thick lines. The crossover to the bare 
gravitational length for large z is confirmed for the high- 
charge particles and was found to be pretty sharp. Con- 
sistently with the theoretical assumption, local charge 
neutrality in the intermediate regime is fulfilled. Second, 
we have tested the location of the transition towards the 
colloidal brazil nut effect by systematically varying the 
mass and charge ratio. The results are summarized in 
Figure |21 The theory, which predicts the transition at 
(Zi -I- 1)7(^2 -I- 1) = TOi/m2, is shown as a straight line 
there. All parameter combinations simulated are indeed 
separated by this theoretical prediction, confirming our 
simple theory. This is remarkable as any wall or bulk 
correlation effects are neglected in the theory. 

In order to elucidate this further, we have compared 
the theoretical predictions 1 /71 , 1 /72 for the heights with 
the simulation data in a situation where the transition 
line was crossed. In Figure the heights are shown as 
a function of a varied charge Zi. The simulation data 
reveal that the variation of both heights hi with Zi is 
large close to the transition and becomes maximal at the 
transition. In the theory, this feature is reproduced and 
accompanied by a generic cusp at the transition. Though 
the cusp is smeared out and in general the heights are 
larger in the simulation data due to the density reduction 
close to the wall, there is still semi-quantitative agree- 
ment. The rapid variation of the heights at the transi- 
tion implies that the location of the transition towards 
the brazil nut effect is very robust explaining the validity 
of the theory in Figure |21 

Finally we address the case of added monovalent salt. 
Density profiles for the two macroion species and the 
counter- and colons are presented in Figure ^ Ad- 
dition of salt reduces the brazil nut effect; in the in- 
set of Figure 01 the heights are plotted versus added 



salt concentration. The threshold salt concentration c. 
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FIG. 4: Density profiles for added salt. The counter- and 
colon densities have been divided by Zi-\-Z2, the colon density 
is smaller than that of the counterlons. The parameters are 
Zi_ = 45, mi/m2 = 2, ^2 = 2Qa, and c^a^ = 2.5 x 10"^. 
Inset: Sedimentation heights hi of heavy (+) and light (0) 
particles as a function of reduced salt concentration for Z\ = 
32, m\/m2 = 2, and ^2 = 5 cr. The theoretical estimate of the 
threshold salt concentration is marked by an arrow. 



needed to reverse the brazil nut effect is estimated as 
c"s^ K, ZiPi{zq) + Z2P2{zo) which is the counterion con- 
centration at the position zq where the two colloidal den- 
sity profiles cross in the salt-free case. This salt concen- 
tration is indicated as an arrow in the inset of Figure 21 
confirming the validity of the estimate. Though the brazil 
nut effect is proven to remain stable with respect to added 



salt, it will only show up for deionized solutions. Highly 
charged suspensions in non-polar solvents with a small 
dielectric constant e and low impurity ion concentrations 
|20j | are promising candidates to exhibit a strong brazil 
nut effect. 

In conclusion, we predict an analog of the granular 
brazil nut effect in equihbrium sediments of charged sus- 
pensions, which is generated by an entropic charge lift- 
ing due to the Coulomb coupling to the counterions. An 
experimental verification of the levitation should be pos- 
sible employing depolarized hght scattering or con- 
focal microscopy |2(J. The simulated charge asymme- 
tries correspond to micelles and proteins rather than 
highly charged colloids but the simple theory which was 
confirmed by the simulations is applicable for arbitrary 
charges. 

The lifting effect is immediately generalizable to sed- 
iments of solutions which are polydisperse in mass and 
charge j^. Therefore the colloidal brazil nut effect has 
important biochemical implications for the separation of 
polydisperse biological matter, such as protein solutions. 
Analytical sedimentation is typically used in an ultracen- 
trifuge to separate different species The colloidal 

brazil nut effect implies that the separation is sensitive 
to the mass per charge but not to the mass itself. 
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